6 Parametric oscillator 



6.1 Mathieu equation 



We now study a different kind of forced pendulum. 

Specifically, imagine subjecting the pivot of a simple frictionless pendulum 
to an alternating vertical motion: 




This is called a "parametric pendulum," because the motion depends on a 
time-dependent parameter. 

Consider the parametric forcing to be a time-dependent gravitational field: 

9{t) =90+ m 

The linearized equation of motion is then (in the undamped case) 

The time-dependence of g{t) makes the equation hard to solve. We know, 
however, that the rest state 

e = e = o 

is a solution. But is the rest state stable? 

We investigate the stability of the rest state for a special case: g{t) periodic 
and sinusoidal: 

9{t) =9o + 9i cos{2ujt) 
Substituting into the equation of motion then gives 

-^+ujl[l + h cos{2ujt)] 6 = (14) 
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where 

ul = go/l and h = qi/qq > 0. 
Equation (14) is called the Mathieu equation. 

The excitation (forcing) term has amplitude h and period 

27r TT 
^ exc ~ 7i ~ 

On the other hand, the natural, unexcited period of the pendulum is 

T 

-'^nat — 



We wish to characterize the stability of the rest state. Our previous methods 
are unapplicable, however, because of the time-dependent parametric forcing. 

We pause, therefore, to consider the theory of linear ODE's with periodic 
coefficients, known as Floquet theory. 



6.2 Elements of Floquet Theory 



Reference: Bender and Orszag, p. 560. 



We consider the general case of a second-order linear ODE with periodic 
coefficients. We seek to determine the conditions for stability. 

We begin with two observations: 



1. If the coefficients are periodic with period T, then if 0{t) is a solution, 
so is e{t-\-T). 

2. Any solution 9(t) is a linear combination of two linearly independent 
solutions 9i(t) and 6*2 (i): 

e{t) = Adiit) + Bd2it) (15) 

where A and B come from initial conditions. (Reason: the system is linear and 
second-order.) 
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Thus 



Since Oi{t) and 62{t) are solutions, periodicity =^ 6i{t-\-T) and 62(1 + T) are 
also solutions. 

Then di{t + T) and 6*2 + T) may themselves be represented as linear com- 
binations of di{t) and d2{t): 

9i{t + T) = aOiit) + (392(1) 
02it + T) = e lit) + 662(1) 

e{t + T) = A[aeiit) + (]e2it)] + B[jeiit) + 6e2it)] 

= {Aa + B'y)ei{t) + {AI3 + B6)e2{t) 
We rewrite the latter expression as 

eit + T) = A'ei{t) + B'e2it) (le) 

where 
or 

a! = Ma, a = ^ ^ ^ • 

Now choose ^ ^ ^ to be an eigenvector of M, with A the associated eigen- 
value. Any other a would have a projection onto one of the eigenvectors. 

Then 

A' = XA and B' = XB. 
Using (16) and (15), we find that 

0{t + T) = XO{t). 

Thus 6{t) is periodic within a scale factor A. The question of stability then 
hinges on the magnitude of A. 
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Define 

(We are interested only in growth or decay, not oscillations due to the exponential multiplier.) 

Then 

e{t + T) = e''^e{t). 

Here 6 is rescaled by e^^ each period. 

Finahy, define P{t) to be periodic such that P{t -\- T) = P{t). The rescahng 
may now be expressed continuously as 

[The above e{t + T) = e'^(*+^)P(t + T) = e''^e''*P(t) = e^''^9{t)] 

StabiUty thus rests on the sign of ji. 

Thus we find that the solution to a linear second-order ODE with periodic 
coefficient — e.g., the Mathieu equation — is of the form 

(exponential growth or decay) * (periodic function of time) 



6.3 Stability of the parametric pendulum 

We proceed to determine under what conditions the rest state of the Mathieu 
equation is unstable, leading to exponentially growing oscillations. 

We expect the tendency toward instability to be strongest when the excitation 
frequency is twice the natural frequency, cjq, of the pendulum. 



This is called parametric resonance. 
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The Mathieu equation (for ^ <C 1) is 



+ cc;5 [l + /icos(2a;t)]^ = 0. 



We take the forcing frequency to be 2a;o + £, £ <C 1: 

e + ul[l + h cos(2u;o + e)t] = 0. 



(17) 



We assume that the excitation amphtude /i <C 1, and seek solutions of the 
general form 

6{t) = a{t) cos ^a;o + t + b{t) sin ^a;o + -e^ t. 

The frequency a;o + is motivated by our resonance argument. Stability will 
depend on whether a{t) and b{t) exponentially grow or decay. 

We begin by substituting 0(t) into the equation of motion. We will then 
retain terms that are linear in e and first-order in h. 

The calculation is messy but is aided by recalling trig identities like* 



cos 



^0 + ^£ j t 



cos [{2cjo + £)t\ = ^ cos 









3 (wo + ^£ j * 


+icos 





higher freq term 



The term with frequency 3((x;o + \e) may be shown to be higher order with 
respect to h. Such higher frequency terms create small amplitude perturba- 
tions of the solution and are neglected. (You shall see on a problem set how 
such terms arise from nonlinearities.) 

Also, we shall retain only those terms that are first-order in e. Thus we 
neglect the 0{e^^ associated with accelerations B. We also assume 

a ~ era, 6 ~ £:6, 

and thereby neglect a, h. These assumptions are validated by the final result. 



*That is, wc use cos A cos B = i[cos(A + + cos(A - S)] and sin A sin B = i[sin(A + + sin(yl - B)]. 
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The result of the substitution is 



— ^2a + b£ + ^huobj uq sin 
+ (2b — ae + ^hujQaj uq cos 



^0 + ^sj t 



= 



For this expression to be true for all t, the coefficients of both cos and sin 
must equal zero, i.e, 

2d + be + l-huob = 
2b — ae + ^hujQa = 



We seek solutions a{t) oc e^^ and b(t) oc e^*. 
Thus, substitute d = /ia, b = jib: 

lia + ]^hu^ = 

- ]^hu^ -^6 = 

These equations have a solution when the determinant of the coefficients of 
a and b vanishes. Thus 



1 



= 



or 



or // = 



(18) 



Instability (parametric resonance) occurs when \i is real and positive. This 
will occur when //^ > 0, or 



43 



This result is summarized in the fohowing phase diagram: 



8 



ho)/2 



exact resonance 





h 



Thus we see that resonance occurs not only for e = 0, but for a range of e 
within the lines ±/ia;o/2. 

The larger the forcing amplitude, the less necessary it is for the force to be 
exactly resonant (cj = 0). 

Conversely, infinitesimal forcing (i.e., h) is sufficient for instability so long as 
6.4 Damping 

The damped parametric pendulum is 



where the factor e^^ results from the periodic forcing. 

The instability boundary is therefore no longer given by = (i.e., equation 



d + 2^d-\-Lul[l-\-h cos((2a;o + e)t)] 6 = 



For unforced oscillations {h = 0), damping produces solutions like 

6{t) ~ e~^* X (oscillation). 

For h > 0, we expect 



e{t) ~ e^^-'^^* X (oscillation) 



(18)), but by 



/i-7 = 0. 
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Instability thus occurs for //^ > 7^. Using equation (18), we obtain 

\ 2 



or 



2 1 

^ =4 



1/2 



< € < 



nl/2 



Setting £ = (exact resonance) shows that instabihty is only possible when 
the quantity in the brackets is positive. Thus we require h > hc^ where 
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/ic = 

The phase diagram is therefore modified accordingly: 




unstable 



h 



N -hco^/l 



6.5 Further physical insight 

We have seen that the instability is strongest at exact resonance. 

That is, if the natural period of the pendulum is Tnat, the most unstable 
excitation has period 

^exc ~ '^'-l^B.at 

corresponding to angular frequencies 

<^exc = 2c<;o 

This is called a subharmonic instability, because instability results at half the 
frequency of the excitation (i.e., ljq = a;exc/2). 
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What other periods of excitation would be expected to lead to instability for 
small h7 

Any integer multiple of Tnat suffices. Why? Because this would correspond 
roughly to, say, pushing a child on a swing every nth time he/she arrives at 
maximum height. 

Therefore (weaker) instabilities will occur for 



Texc = integer x 
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